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Different criteria of constructing clusters and tracing back A resonances from the intermediate- 
energy neutron-rich HICs are discussed by employing the updated UrQMD transport model. It is 
found that both the phase-space and the coordinate-density criteria affect the single and the double 
neutron/proton ratios of free nucleons at small transverse momenta, but the influence becomes 
invisible at large transverse momenta. The effect of different methods of reconstructing freeze-out 
1 As on the A°/A ++ ratio is strong in a large kinetic energy region. 
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^ 1 I. INTRODUCTION 

< 

I/"") \ The equation of state (EoS) of nuclear matter is one of the most important topics in the field of low and intermediate 
energy heavy ion collisions (HICs). In recent years, more and more work have been focusing on the isospin asymmetry 
of the HICs, which introduces more uncertainties into the EoS. Without considering the momentum dependence of the 
EoS, the isospin dependent EoS for asymmetric nuclear matter can be simply expressed as e(u,S) = eo(«) + e sym S 
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where u = p/po is the reduced nuclear density and 5 — (p n — p p )/p is the isospin asymmetry in terms of neutron and 
proton densities, eo and e sym are the isospin independent term (which includes the Skyrme and the Yukawa potentials) 
and the symmetry energy (which can be expressed as e sym = SqF(u) where So is the symmetry energy coefficient 
and F(u) the density dependence), respectively (see, e.g., Refs. [Hll]). Since the symmetry energy term plays also an 
important role in nuclear structure and astrophysics, it indeed deserves more attention. With the many efforts from 
both theoretical and experimental sides, the uncertainties in EoS have been largely constrained although arguments 
remain in this field. For example, (1) the isospin-independent EoS has been constrained into a soft compress modulus 
' (K ~ 190 — 270 MeV, see, for example, [HUH, Q)- (2) The symmetry energy coefficient has been constrained into 
. the region So — 30 — 36 MeV 0, 0, @|- And, (3) the strength factor (7) of the density dependence of the symmetry 
potential is also shown not to be stiff at subnormal densities i.e., 7 = 0.69 or 1.05 when the form F(u) = u 1 

is employed, which depends on the treatments of nucleon-nucleon collisions in the BUU transport model. 

While the density dependence of the symmetry potential at supranormal densities is still quite uncertain. So 
far several (probably) sensitive probes especially on this issue have been brought out by theoretical groups, which 
are hoped to be observed from experiments such as the planned FRIB(USA), the upcoming GSI new facility FAIR 
(Germany), and the cooling storage ring CSR being tested at Lanzhou (China), etc. For example, the w~ to 7r + 
• i-H , multiplicity ratio and the neutron-proton differential collective transverse flow were firstly pointed out in Ref . . 
While the threshold production of pions and kaons and the difference between neutron and proton elliptic flows were 
stressed in Refs. [H,|r|. In our previous investigations, we found that the £~/E + ratio [3] at threshold, the n + — n~ 
ellip tic flow difference at moderate transverse momentum |15| . and the A~/A ++ ratios at large transverse momenta 
[la ] are sensitive to the density dependence of the symmetry potential at high densities. Despite of the difficulties of 
detecting and analyzing A resonances, the advantage of taking A as a probe for symmetry potential at high densities 
is also obvious: firstly, it is known that most of As can be produced from the high density region. Secondly, the 
A t> Nir loop sustains the A matter for a relatively long time, hence the evolution of As should be heavily influenced 
by the mean field. Finally, since some of nucleon- and the pion-related quantities are predicted to be sensitive to the 
symmetry potential, it is believed that the corresponding A-related quantities will show strong effect on symmetry 
potential as well. Thus we are interested in checking one crucial question: can the transverse momentum distribution 
of the A~/A + ratio, or more realistically, the kinetic energy distribution of the A°/A + ratio (since it is not easy 
to reconstruct A - and its momentum components through its neutron and n~ daughters), be experimentally taken 
as a sensitive candidate for detecting the symmetry potential at supranormal densities? Besides, we found that the 
transverse momentum dependence of the neutron/proton (n/p) ratio of free nucleons is also sensitive to symmetry 



* e-mail address: liqf@fias.uni-frankfurt.de 



2 



potential fl7j | . However, since the effect of symmetry energy term is secondary when comparing it with the isospin- 
independent term, one should pay attention to the late stage of HICs in any transport-model related calculations. 
In this work, we also intend to test the influence of different criteria of constructing clusters on the n/p ratio of free 
nucleons. 

The paper is arranged as follows: in section 2, the UrQMD model is introduced briefly and the projectile-target 
combinations to be investigated are chosen. Section 3 explains different criteria of constructing free nucleons/fragments 
and tracing back A resonances at freeze-out. The effect of the reconstruction criteria on the multiplicities of nucleons 
and As is discussed. In section 4, the influence of different reconstruction criteria on the transverse momentum 
dependence of the n/p ratio and the kinetic energy dependence of the A°/A ++ ratio is illustrated and discussed. The 
conclusions are given in section 5. 

II. SYSTEM DETERMINATION 

The UrQMD transport model [H, EH has been updated for the investigations of intermediate energy HICs. For 
details of the updated version of the UrQMD model, the reader is referred to Ref. |20| and references therein. In 
this work, the soft EoS with a momentum dependent term (SM-EoS) is adopted. Two types of density dependent 
symmetry potential energy are selected: (1) F(u) — u 1 , where 7 factor varies from 0.5 (dubbed as "g05", soft), 
1.0 ("gl", moderate) to 1.5 ("gl5", hard). (2) the so-called "DDH/?*" (very soft) symmetry potential energy from 
the relativistic mean- field calculation [2l| . Furthermore, The "Dirac"-type medium modification of nucleon-nucleon 
elastic cross sections discussed in [22j is also considered here. 

The choice of the reaction system should be with great care as well. For A-related probes, first of all, in order to 
excite enough A(1232) resonances but not to excite too many other higher-lying resonances, the HICs at moderate 
beam energies such as 400 ~ 1000 A MeV should be chosen. At the same time, considering of the large background 
spectrum of non-correlated (p,7r) pairs which is induced especially by heavy system collision (23|, a midsize and 
neutron-rich projectile-target system is necessary in order also to obtain an initially large isospin asymmetry S value. 
Further, inspired by the experiments with 2 or 4 sets of projectile-target systems of different asymmetry [24], [25| and 
the model-based investigations [2(1 Hfij] in order to cancel other potential effects which are not (obviously) related 
to isospin, we select the central (6 = 0-2 fm) 132 Sn + 132 Sn and 112 Sn + 112 Sn reactions at E b = 800A MeV as 
examples. While for nucleon-related probes, the HICs (besides the reactions with respect to Sn-isotopes we select the 
Pb+Pb system as well) at a lower beam energy such as 400A MeV are preferred and adopted in this work. 

III. DIFFERENT CRITERIA OF CONSTRUCTING CLUSTERS AND TRACING BACK A 

RESONANCES 

Generally speaking, two methods for constructing clusters exist in the transport model calculations. The one 
is called as the coalescence model, which is normally used in the QMD-like analysis. In the coalescence model 
nucleons with relative distances |Ar| < Ro and relative momenta |Ap| < Po are considered to belong to one cluster 
(depicted in the left plot of Fig. [1] titled as "T-I"), otherwise, the nucleons are free. The other one adopts a 
coordinate density cut p c (depicted in the right plot of Fig. [TJ titled as "T-II" ) to separate " free" nucleons from 
other fragments. This method is often used in the BUU-like model analysis, in which the density of each real particle 
is obtained with respect to all of its test particles. While in QMD-like transport models, each particle is represent 
by a Gaussian wave packet in the phase space, the density of each particle in coordinate space can be expressed as 

< r i- r j) 2 

Pi = J p(ri)pdr = J p(r\) ^V p(rj)dr = ^^3/2 X)j e 51 ■ I n this work, we examine the freeze-out condition 
Pi < Pc = Po/10 (dubbed as "Rc") to "find" free nucleons. Fig. shows the time evolution of the fraction of nucleons 
to be free from Pb+Pb central collisions at Et, = 0.4A GeV (by averaging 100 events). The SM-EoS with a linear 
symmetry potential ("gl") is adopted as an example. In the "T-I" mode, the free nucleons are found out with two 
sets of R and P parameters: (1) "RP-I": R = 2.8 fm and P = 0.2 GeV/c, (2) "RP-II": R = 3.5 fm and P = 0.3 
GeV/c. In the "T-II" mode, the "Rc" parameter set is used. It is easy to understand that with the increase of Ro and 
Po values, the number of free nucleons becomes less and less. Meanwhile, it is seen that the number of free nucleons 
with "RP-I" mode is almost same as that with "Rc" mode after t ~ 80fm/c and not sensitive to the time evolution 
any more. Before this time, much more nucleons in "RP-I" mode than in "Rc" mode are taken as free due to the 
large momentum difference but the relatively small distance difference between nucleons. Therefore, the choice of the 
cut-time t c to stop the transport program will affect the final freeze-out of nucleons if a too short time is selected. 
From now on, the t c = lOOfm/c is selected for further investigations. 

In order to track the resonances at freeze-out, one needs to look through the detailed history of the time evolution 
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FIG. 1: Two methods for constructing clusters. Left plot: the relative distance of two particles |Ar| < Ro and the relative 
momentum of two particles |Ap| < Po ("T-I"), Ro and Po are free parameters in units of fm and GeV/c. Right plot: coordinate 
density cut for the particle i in fragments: pi > p c — po/c ("T-II"). 
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FIG. 2: Time evolution of the fraction of "free" nucleons from Pb+Pb central collisions at Et = 0.4A GeV. The SM-EoS with 
a linear symmetry potential is adopted in calculations. No other physical cuts are selected. Results of "finding" free nucleons 
from different criteria ("RP-I", "RP-II", and "Rc") are compared (see context). 
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FIG. 3: Two methods for reconstructing A resonance at freeze-out. Left plot: the pion produced from a A decay is found and 
the pion does not rescatter any more. It is not concerned if the other nucleon daughter collides further with other particles 
or not (marked with a shadowed box). This mode is called as "M-I". Right plot: only the As whose decay daughters do not 
rescatter any more are selected. This mode is called as "M-II" . 
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FIG. 4: Kinetic energy distribution of A and A ++ resonances at freeze-out which is "found" by the two methods depicted in 
Fig.rjfrom central 132 S7i+ 132 Sn collisions at E b = 0.8A GeV. The SM-EoS with a stiff symmetry potential ("gl5") is adopted. 
The rapidity cut of pions < 0.5 is used. 



5 



of HICs. In the UrQMD model there is a standard OSCAR-formatted output file ("ftn20") which includes complete 
event history, that is, besides the initial and the final freeze-out states, all binary collisions, string- fragmentations, 
and hadronic decays of particles are recorded into this file. From this file, one can easily find the A decay channel 
A — > Ntt. It is known that the produced daughters nucleon and pion from this channel will probably collide further 
with other particles [27| . Thus it is not easy to detect A resonances at freeze-out time experimentally. Currently there 
exist two methods in experiments which are employed to "find" baryon resonances: (I) the one-pion pt distribution 
is defolded to yield the mass distribution; (2) the correlated proton and charged pion pairs are analyzed to yield the 
invariant mass distribution [23]. In Fig. [3] we also show two methods for "finding" A resonance at freeze-out: (I) in 
"M-I" , the A in which decay channel the pion daughter will not further rescatter is found out while, it is not concerned 
if the other nucleon daughter rescatters or not. This mode can be taken as a "contaminated" reconstructing mode. 
(2) in "M-II" , both the pion and the nucleon from the decay of the A resonance at freeze-out are confirmed, the A 
resonance is thus taken to be reconstructible. And this mode is a "clean" reconstructing mode. Fig. 2] illustrates the 
kinetic energy E^in distribution of A and A ++ resonances at freeze-out "detected" by the above two methods for 
central 132 Sn + 132 Sn collisions at E b = 0.8A GeV. The SM-EoS with a stiff symmetry potential ("gI5") is adopted 
as an example. The rapidity cut of pions \y n \ < 0.5 (y = \ \og( ^ 1Tl+p " ), E cm andp// are the energy and longitudinal 

momentum of the pion in the center-of-mass system) is used. With the "M-I" method the amount of As at small 
kinetic energies is much larger than that with the "M-II" method, which is due to the high chance of rescattering of 
the produced nucleon with others. The difference of the A yields with the two methods is seen to disappear at large 
Ekin due to the fact that the rescattering probability of the fast nucleons becomes rather rare. 



IV. THE INFLUENCE OF DIFFERENT RECONSTRUCTION CRITERIA ON SENSITIVE 

OBSERVABLES 

Fig.[5]shows the transverse momentum pt dependence of the n/p ratio of free nucleons from central Pb+Pb collisions 
at Eb = 0.4A GeV. Two sets of density dependent symmetry potential "gl5" (hard) and "DDHp*" (very soft) are 
adopted in calculations. And three modes "RP-I", "RP-II", and "Rc" are taken into account for constructing clusters 
at the cut-time t c = fOOfm/c. No rapidity cut is selected. First of all, it is known that the n/p ratio at low and high 
transverse momentum reflects the behavior of symmetry potential at subnormal and supranormal densities, which has 
been systematically studied in our previous work [17|, in which the "RP-II" is used in the coalescence-model analysis. 
With shorter relative phase-space distances used in the "RP-I" mode, more nucleons are counted into free nucleons, 
which drives the n/p ratio at small pt to approaching the initial n/p ratio of the system (=1.54). Surprisingly, although 
the total amount of free nucleons from "RP-I" mode is seen almost same as that from "Rc" mode, the n/p ratio of 
free nucleons is visibly different from each other at small pt- the "Rc" mode gives a larger n/p ratio. It implies that 
the two methods of constructing clusters used for transport models are not equivalent at small pt and bring more 
uncertainties when one hopes to evaluate the sensitivity of nucleon-related probes to the symmetry potential. One also 
finds that the softer the symmetry potential is, the more sensitive to the free n/p ratio the methods of constructing 
clusters are. It is easy to understand since the soft symmetry potential has larger value than the stiff one at subnormal 
densities, thus the softer symmetry potential influences the freeze-out of particles more obviously at the late stage. 
Then, let us move on to see the n/p ratio in the large pt region (p t > 0.6GeV/c), it is interesting to see that all of 
the three constructing modes have almost no different effect on the values of the n/p ratio. It is obviously due to the 
fact that the free nucleons with large pt are dominantly emitted from high densities at early stage and are influenced 
very weakly by the late stage. It also shows that the n/p ratio of free nucleons at large pt is a (relatively) clean and 
sensitive probe to the symmetry potential at high densities. 

The pt dependence of the n/p ratio from 132 Sn + 132 Sn and 112 Sn + 112 Sn systems are also calculated. The results 
from the two systems are shown in the left plot of Fig. [5] separately. In this figure, the results with modes "RP-I" and 
"Rc" are compared and a rapidity cut \y\ < 0.2 is employed. The results are similar to those in Fig. [5] The effects of 
the symmetry potential and the freeze-out criteria on the n/p ratio are week in the more isospin-symmetric system 
112 Sn + 112 Sn. The right plot of Fig. [5] illustrates the double ratio expressed as the n/p ratio from 132 Sn + 132 Sn 
system divided by the n/p ratio from 112 Sn + 112 Sn system, and is dubbed as {n/p)sni32/(n/p)snii2- Similar to the 
single n/p ratio shown in the left plot of Fig.[6j the effect of freeze-out criteria is still seen at small pt, especially for 
the soft symmetry potential case, while at pt > 0.6GeV/c, it almost disappears. Therefore, both the single and the 
double ratios at large transverse momenta are visibly affected by the density dependence of the symmetry potential 
but invisibly influenced by the different criteria of constructing clusters and can be taken as good observables for 
symmetry potential at high densities. 

Let us finally check the effect of the different analyzing methods of reconstructing A resonances at freeze-out on its 
particle ratio A°/A ++ . In Fig. [7] (left plot) we present the kinetic energy (Ef-in) distribution of the A°/A ++ ratio at 
freeze-out where As are collected with the methods "M-I" and "M-II" introduced above. The symmetry potentials 
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FIG. 5: Transverse momentum p t dependence of the n/p ratio for central Pb+Pb collisions at Eh = 0.4A GeV. Two sets of 
symmetry potential "gl5" and "DDHp*" are adopted in calculations. Three modes "RP-I", "RP-II", and "Rc" are taken into 
account for constructing clusters at the cut-time t c — lOOfm/c. 
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FIG. 6: Left plot: Transverse momentum p t dependence of the n/p ratio of free nucleons for central 132 Sn+ 132 Sri and 
112 Sn+ 112 S'n collisions at = 0.4A GeV. Two sets of symmetry potential "gl5" and "DDHp*" are adopted in calcula- 
tions. Two modes "RP-I" and "Rc" are taken into account to construct clusters. The rapidity cut \y\ < 0.2 is chosen. Right 
plot: The double ratio between n/p ratios from the 132 Sn+ 132 5n system and from the 112 Sn+ 112 5'n system. 
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FIG. 7: Left plot: Kinetic energy distribution of the A°/A ++ ratio at freeze-out. The As are determined with two methods 
"M-I" and "M-II" (see context). The central 132 Sn+ 132 S'n collisions at E b = 0.8A GeV are calculated. The rapidity cut of 
pions \y„\ < 0.5 is chosen. Right plot: The A°/A ++ ratio from 132 Sn+ 132 SVi is scaled by the A°/A ++ ratio from 112 Sn+ 112 SVi. 



"g05", "gl5", and "DDHp*" are adopted. For each case 0.36 million central 132 Sn+ 132 Sn collisions at E b = 0.8A GeV 
are calculated. And a rapidity cut of pions \y n \ < 0.5 is adopted as well. First of all, the "crossing" behavior of the 
ratios with soft and stiff symmetry potentials can hardly be seen from the "M-I" reconstruction mode. It implies that 
the A°/A++ ratio versus E^m is sensitive to the density dependent symmetry potential at high densities. However, 
in the "clean" "M-II" mode, the "crossing" behavior re-occurs which means the symmetry potential at low densities 
shows its role. Actually, the two modes give a quite different A°/A ++ ratio especially at the small Ekin- Even at 
the large Eki n (> 0.5GeV), the difference does not disappear. The "clean" mode "M-II", which implies that the 
selected As come from more outer space where densities are low, gives a much large ratio at all Ekin due to the fact 
that these As are produced from more neutron-rich region. Even in the scaled ratio by the ratio from 1 12 Sn-I- 1 12 Sn 
system (dubbed as (A°/A ++ )s„i32 — (A°/A ++ )s„ii2) shown in the right plot of Fig. [7| the huge effect of different 
reconstructing criteria on the scaled ratio still exists. Although the A°/A ++ ratio indeed shows its sensitivity to the 
density dependence of symmetry potential, the detailed process of the reconstruction for As at freeze-out should be 
paid much more attention. 

In summary, different criteria of constructing free nucleons and tracing back A resonances from the intermediate- 
energy neutron-rich HICs are discussed with the help of the UrQMD model. It is found that both the parametrization 
of relative phase-space distances used in the coalescence model and the coordinate density cut for free nucleons modify 
visibly the single and the double n/p ratios at small transverse momenta, but the influence is weak at large transverse 
momenta. Hence, the single and the double n/p ratios of free nucleons at large transverse momenta can be taken 
as sensitive probes for symmetry potential at supranormal densities. The momentum dependence of the symmetry 
potential might complicate this situation and deserves further careful investigation [28| . Methods of reconstructing A 
resonances affect the A°/A ++ ratio versus kinetic energy strongly, which should be paid more attention if one wishes 
to employ it as a probe to detect the symmetry potential of the nuclear matter at high densities. 
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